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Scaling up linear-optics quantum computing will require multi-photon gates which are compact, 
phase-stable, exhibit excellent quantum interference, and have success heralded by the detection of 
ancillary photons. We investigate implementation of the optimal known gate design which meets 
these requirements: the Knill controlled-Z gate, implemented in integrated laser-written waveguide 
arrays. We show that device performance is more sensitive to the small deviations in the coupler 
reflectivity, arising due to the tolerance values of the fabrication method, than phase variations in 
the circuit. The mode fidelity was also shown to be less sensitive to reflectivity and phase errors 
than process fidelity. Our best device achieves a fidelity of 0.931 ±0.001 with the ideal 4x4 unitary 
circuit and a process fidelity of 0.680±0.005 with the ideal computational-basis process. 


I. INTRODUCTION 

Recent advances in on-chip integration of efficient pho¬ 
ton sources [1-4] and detectors [5-7] hold promise for 
significant progress in the linear optics architecture for 
quantum information processing and simulation. This 
architecture relies on the insight that entangling quan¬ 
tum gates can be realised probabilistically by interacting 
photonic qubits using effective optical nonlinearities in¬ 
duced by measurement [8]. However, efficient scaling in 
this architecture necessitates entangling gates which are 
both logically and physically scalable. Logical scalabil¬ 
ity requires successful operation to be heralded nonde- 
structively, typically by the detection of additional ‘an- 
cilla’ photons. Physical scalability, meanwhile, requires a 
compact and phase-stable architecture. Integrated, non- 
heraided entangling gates have been demonstrated [9, 10] 
as have heralded gates in bulk optics [11, 12]. However, 
gates which are scalable in both the logical and physical 
senses have not been implemented to date due to their 
geometric complexity as well as the need for low circuit 
losses. 

The simplest heralded entangling two-qubit photonic 
gate design, with the highest known success probability, 
was found by Knill [13] and implements a controlled-Z 
operation with probability 2/27. This heralded CZ design. 
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henceforth called the hcz, relies on pairwise non-classical 
interference of four indistinguishable photons in a cir¬ 
cuit with four particular beamsplitters (BSs), as shown 
in Fig. 1(a), as well as a stable phase shift of precisely tt 
between the first and second BS pairs. 

Integrated arrays of coupled waveguides could enable 
compact, phase-stable circuits with the requisite split¬ 
tings and phase for Knill’s design. However, integrated 
hcz circuits require a physical swapping of neighboring 
modes; such crossovers are difficult to achieve lithograph¬ 
ically. Laser-written waveguides have recently been used 
to demonstrate a wide variety of quantum photonics cir¬ 
cuitry [14, 15]. The femtosecond-laser direct-write tech¬ 
nique in particular allows 3D waveguide arrays, simplify¬ 
ing waveguide crossovers, along with demonstrated high 
mode indistinguishability [14], and has recently been em¬ 
ployed for multiport and arbitrary-phase directional cou¬ 
plers [16-18], all-optical routers [19], circuits for small- 
scale quantum simulations [20-22], quantum walks [23], 
and non-heralded quantum gates [10]. 

Here we investigate implementation of the hcz gate 
both theoretically and experimentally using the direct- 
write technique, with a particular focus on its ac¬ 
tion in the presence of deviations from optimal phase 
and reflectivity parameters. We derive the variation 
in gate performance as quantified by two metrics—the 
optical-circuit mode-fidelity and the computational-basis 
process-fidelity—with respect to such deviations. We fur¬ 
ther detail the fabrication of 12 prototype circuits in¬ 
cluding a novel and simple method for achieving the req¬ 
uisite internal phase, and their full characterisation us- 
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ing coherent techniques as well as quantum interference 
which confirms their excellent mode indistinguishability 
and suitability for the single-photon regime. 

II. MATERIALS AND METHODS 

The optical circuit for the hcz gate is shown in 
Fig. 1(a). The control and target qubits are each en¬ 
coded as single photons across a pair of modes, labelled 
C and Co for the control, and T and To for the tar¬ 
get. Co and To, which encode |0) for their respective 
qubits, interact with neither the |l)-modes nor the an- 
cillas; our fabricated circuits contain only the four in¬ 
teracting modes as shown in Fig. 1(b). Modulo local 
phases on the input and output modes, this circuit im¬ 
plements a heralded CZ operation: conditioned on the 
detection of one photon in each ancilla mode it flips the 
sign of the |ll)-term of an arbitrary two-qubit input state 
0^00 | 00 ) + 0^10 | 10 ) + 0^01 | 01 ) + 0^11 | 11 ). 

A. Circuit modelling and design. 

Demonstrations of entangling linear optics quantum 
gates relying on post-selection of the computational pho¬ 


tons go back more than a decade and include both free- 
space [24, 25] and integrated implementations [9, 10]. 
In such gates, success relies on detecting the photonic 
qubits in particular output modes, precluding their use in 
subsequent multi-qubit operations [ 8 ], whereas heralded 
gates—where success is signalled by detection of ancil¬ 
lary photons without measuring the output qubits—can 
be incorporated as modules in complex quantum com¬ 
putations. The effects of fabrication imperfections on 
post-selected gates have been modelled [26, 27]. How¬ 
ever such gates can effectively be simplified to include 
only a single instance of two-photon interference at one 
beamsplitter [25]; by contrast the the /iCZ, despite its 
simplicity relative to other heralded gates, requires four 
distinct two-photon interference events at four beamsplit¬ 
ters. Here we quantify the effects of device imperfections 
on such a complex heralded gate. 

We first modei imperfect hcz device operation as foi- 
iows. We assume four single-mode waveguides coupled 
by BSs as in Fig. 1(b). The quantum state of the light 
is described by four bosonic creation operators a^, a^, 
a^, and which create a photon in the control, target, 
and two ancillary modes respectively. We employ the 
symmetric BS convention so that two modes ai and 02 
transform as 


BS((9): 


( 1 ) 


We first consider a photon or coherent state in a super¬ 
position of the input modes into a circuit of the form in 
Fig. 1(b). Our modei neglects intrinsic loss, as in our 
fabricated circuits such loss is nearly constant across the 
waveguides and therefore results only in a reduction in 
the output state amplitude and, thus, the gate success 
probability. The circuit then maps the input creation 
operators to the outputs via the transformation 
where is a unitary matrix, and the 
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mode indices k and j are ordered {C, T, A, T}. We allow 
for arbitrary splitting parameter angles On^ ^={1, •••5 4}. 
The requisite internal phase shift is implemented by an 
additional phase of tt on BSS, equivalent to Os —Os- 
Using Eq.(l) and allowing for additional undesired phase 
shifts (aj^ ^ between the BS pairs, we find that 

all these unwanted internal phases can be collected into 
a single net phase shift For the to¬ 

tal circuit action, modulo external local phases, we thus 
find: 


cos 0 i cos Os cos O 2 sin Os cos 0 i sin Os sin O 2 sin Os 
cos 0 i sin Os — cos O 2 cos Os sin 0 i sin Os — cos Os sin O 2 
cos 6^4 sin 6 >i sin 6^2 sin 6^4 — cos 6 >i cos 6^4 cos ^2 sin 6^4 

sin 0 i sin 0 ^ cos 6^4 sin O 2 — cos 0 i sin O 4 cos O 2 cos 6^4 


Up to external phases and in the absence of net phase 
(j)N^ the ideal matrix given by Knill [13] is achieved 
by the target angles of 0 i= 02 = 0 s= arccos \/l/3, and 
6 > 4 = arccos 

We employ two metrics to assess the design accuracy 


of a physical circuit for a hcz gate. The mode fidelity 
Fm directly compares the 4x4 circuit mapping matrix 
to the ideal unitary matrix and is given 

by the normalized Hilbert-Schmidt inner product [28]: 
Fm=|Tr{U^''^t[/circ|| 2 /^ 2 ^ A = 4 is the number 

of modes. This metric most closely captures the differ- 
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FIG. 1. (a) The circuit for a hcz gate showing paths for ancillary photons A and B as well as the computational qubits; the 
control (target) photon is encoded across spatial paths Co (To) representing |0) and C (T) representing |1). The |0)-modes do 
not interact in the gate; the four remaining modes undergo four beamsplitting operations with reflectivities R(0n)= cos^(On) 
as described in Eq. (1). The light-coloured side indicates the surface yielding a relative tt phase change upon reflection, (b) 
The four interacting circuit modes modelled as a waveguide array, showing the crossover and optimal reflectivities for the BSs 
implemented using evanescent coupling. The waveguides are separated by 127 /xm at the device end facets; fan-in and fan-out 
regions are not shown. The input mode labeling is reversed compared to (a) due to the reflectance of the couplers being defined 
as the proportion of input light which couples from one waveguide to the other. 


ences between the manufactured integrated device and 
the ideal target device, but only partially captures how 
the device would function with qubits since it ignores the 
effect of measurement heralding. For instance if the an- 
cilla modes were swapped prior to detection, Fm would 
decrease without any operational effect on the function of 
the heralded gate on the qubits, as it is irrelevant which 
detector detects which ancilla. The second metric di¬ 
rectly assesses the effect of the measurement-induced op¬ 
tical nonlinearity in the space of the qubits, and thus will 
be independent of such irrelevant changes. A quantum 
process S can be represented abstractly as a quantum 
state ps via the Jamiolkowski isomorphism [29], and the 
natural figure of merit for gate quality is then the pro¬ 
cess fidelity: Fp=Ti{p\.r^ps} which simply compares the 
state representing the implemented process and the state 
Pcz representing the process for an ideal CZ gate [30, 31]. 
While we would expect the two measures to be roughly 
similar, particularly for small imperfections in e.g. split¬ 
ting ratios, there is in general no simple relationship be¬ 
tween them, as they are differently sensitive to imperfec¬ 
tions. 

The most direct way to calculate pg is to con¬ 
sider a maximally-entangled state |0max) between the 
Hilbert space on which the process acts, and an¬ 
other fictitious space of the same dimension. The 
process acts on one half of the entangled state, and 
the resulting total state is exactly pg. As our com¬ 
putational input is two qubits, the entangled state 
is |0max) =(|0000) + 10101) + |1010) + |llll))/2. After 
a CZ operation on the first two qubits the result 
is |(/)cz)=(|0000) + |0101) + |1010)-|llll))/2, and the 
corresponding final state representing the process is 
pcz— |0cz) (0CZ I • 

Given that each qubit comprises a photon in two 
modes, |0max) involves four photons encoded across eight 


modes, where the fictitious additional control (target) 
mode has creation operator a ^2 (^t 2 )' With the ad¬ 
dition of the two ancillary modes, the entangled input 
state is thus represented using boson creation operators 

as (l+a^a^ 2 +<^c'^c' 2 +^c^C' 2 ^T^T 2 )^A^B 1^) where | 0 ) 
is a multimode bosonic vacuum, and creation operators 
for the non-interacting logical |0) modes are again omit¬ 
ted. The circuit transforms a^, a^, a^, and accord¬ 
ing to and gate success is heralded by measuring 

exactly one photon in each ancillary mode. This mea¬ 
surement removes these modes and induces a CZ on the 
remaining modes. 


A subtle problem arises when the photonic gate is not 
perfectly balanced. There is then a non-zero amplitude 
for the states proportional to (a^)^ and (a^)^, which lie 
outside the qubit space and represent errors. In charac¬ 
terising circuit performance, we account for these errors 
by calculating the process fidelity against a version of pcz 
which is extended to include these two states but with 
zero support, and thus any weight on these terms will 
always reduce Fp. 


Figure 2 shows the variation of the process and mode 
fidelity due to a deviation in a single BS reflectivity or 
the net internal phase. Both fidelity metrics are much 
less sensitive to small deviations in phase than in split¬ 
ting ratios. Fig. 2 also shows two fidelity distributions 
each resulting from 2000 randomly-chosen gate simula¬ 
tions with simultaneous deviations in all five reflectivity 
and phase parameters. Perhaps unsurprisingly the mode 
fidelity is far less sensitive to errors overall. 
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FIG. 2. Variation of the model mode fidelity Fm and process fidelity Fp with deviations A from the ideal BS angles and internal 
phase shift. The ideal phase is zero and A^n represents any net extra phase introduced between beamsplitters. For the BSs 
A is the variation from the ideal angle; the total reflectivity will be cos^(^ideai + AO). In both cases A is a length variation in 
the physical device. The top graphs show the hdelity when one A parameter is varied and the rest are held at zero. The points 
shown represent the deviations found in our best experimentally-characterised circuit; see Fig. 3 for further details. Note that 
on this scale both the curves and points for BSl and BS2 are indistinguishable. The bottom graphs show the hdelities for 2000 
simulated instances of hcz gates with all A parameters drawn randomly from Gaussian distributions with mean 0 and standard 
deviation a a- The purple distribution has a a = 0.1, which is similar in magnitude to most of the A parameters from our best 
measured circuit; the green distribution has a a = 0.05 in order to show the fidelities achievable with a modest improvement in 
fabrication accuracy. The resulting green (purple) distributions have means of 0.994 (0.980) for Fm and 0.962 (0.863) for Fp. 
Note the logarithmic scale on the horizontal axes. 


B. Device fabrication 


The circuits were fabricated using the FLDW tech¬ 
nique wherein a tightly focused femtosecond laser gen¬ 
erates a localised refractive index contrast in a glass sub¬ 
strate. By translating the glass in {x, y, z) with respect 
to the incident laser, arbitrary 3D regions of net-positive 
refractive index change can be produced. Our fabrication 
employed a titanium sapphire oscillator (800 nm centre 
wavelength, < 50 fs pulse duration) with 5.1 MHz rep¬ 
etition rate [32, 33]. A telescope was used to overfill 
the input pupil of a 100 x oil immersion objective which 
focused the laser into the boro-aluminosilicate sample 
(Corning Eagle 2000) for writing with 66 nJ pulses at 
a translation speed of 1200 mm/min. The sample was 
subsequently annealed to obtain a more symmetric and 
Gaussian refractive index profile [34]. This significantly 
improves waveguide throughput efficiency, as shown in 
(ref. 35). This process yields a mode field diameter of 
5 /im which has excellent overlap with an 800 nm single¬ 
mode optical fibre, and gives rise to fibre-to-fibre cou¬ 
pling loss of only 1.8 dB for straight waveguides of length 
40 mm. Input- and output-coupling accounts for the ma¬ 
jority of this loss, while intrinsic propagation loss is be¬ 


low 0.2 dB/cm [35], and importantly intrinsic losses were 
found to be constant across the four waveguides to within 
measurement error. 

The splitting ratios of the waveguide BSs—or di¬ 
rectional couplers—can be adjusted by changing their 
coupling lengths. Using a symmetric phase conven¬ 
tion, the amplitudes in an ideal directional coupler 
of total length L with uniform coupling constant C 
vary sinusoidally with propagation length ^ as aj ^ 

at cos(Cz)+ia^ sin(C 2 ;), where j,/c={l, 2 } and 0 < 2 ; < 
L. Since the waveguides have nominally identical pro¬ 
files, the reflectivity takes the simple form R = cos^ 7 , 
where 7 = C{z)dz and we allow for variation in the 
coupling strength C{z) along the waveguide [36]. 

While specific, arbitrary phase shifts are difficult to re¬ 
alise precisely without active elements using FLDW [18, 
22 ], adjustments in coupler length also allowed us to 
achieve the requisite internal phase shift of tt. Extend¬ 
ing L such that 7 goes from 0 to {27r—0) changes the 
action of the splitter to aj ^ aj cos(—^)+m^ sin(—^) = 

at cos( 6 >) — m^sin( 6 >). Exploiting this identity, we im¬ 
plemented the requisite phase shift by lengthening BS3 
from 7 = arccos y/l/3 to 7=(27r— arccos y/l/3). Eor ideal 
couplers the relative phase is limited to ±7r/2 and the 
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application of this technique on BS3 yields no undesired 
internal phase even for slight errors in L. In prac¬ 
tice 0Ar 7 ^ 0 can occur due to slight variations in local 
waveguide profile resulting from laser power fluctuations 
in fabrication, as well as from small internal path length 
variations. 

An extensive parameter study of directional couplers 
was completed in order to determine the optimal laser 
characteristics, writing algorithm, and coupling lengths 
for achieving the desired reflectivities and internal phase. 
However, the performance of couplers written according 
to a particular algorithm will nevertheless vary from sam¬ 
ple to sample, depending on the precise substrate and 
laser characteristics at the time of fabrication. In partic¬ 
ular, slight refractive index differences between the two 
waveguides constitute a significant source of deviations 
from intended reflectivities. Such differences yield phase 
mismatch which prevents full power transfer between the 
waveguides, an effect which becomes more pronounced 
as coupling length increases. Twelve separate candidate 
circuits were thus fabricated, both to increase the like¬ 
lihood of achieving near-optimal phase and reflectivity 
parameters in one or more circuits, and in order to ex¬ 
perimentally investigate the sensitivity of device opera¬ 
tion to parameter variations. The origin of the parameter 
variations can be traced to the fact that the writing laser 
is passively mode locked and the cavity is very long in or¬ 
der to achieve a balance of sufficiently high pulse energy 
(>100 nJ) combined with high repetition rates (5 MHz). 
This cavity is 30 m long and any perturbation can have 
a substantial knock on effect on the resultant pulse en¬ 
ergy at short (sub second) time scales. A change in the 
pulse energy by as little as 5% can result in a change 
in refractive index variation and directional couplers are 
extremely sensitive to these variations. In principle these 
laser fluctuations could be improved by temperature sta¬ 
bilisation of the cavity area or potentially operation of 
the cavity in-vacuuo. 


III. RESULTS AND DISCUSSION 
A. Coherent device characterisation. 

We characterised the fabricated candidate circuits us¬ 
ing a recently-demonstrated technique [37, 38] which 
yields = r™®^®exp(i(/)™®^®) using only single- and 

two-mode bright coherent states and output intensity 
measurements. The moduli result from intensity 

measurement at each output k for an input at mode j. 
The phases are obtained as follows: a two-mode 

coherent state is injected into two inputs, and a relative 
phase between the modes is induced via continuous path- 
length variation in one mode using motorised translation. 
The output interference fringes are recorded with fast 
photodiodes and an oscilloscope, and the phases 
are simply the phase differences between the pairs of re¬ 
sulting periodic output intensity signals How¬ 


ever, due to experimental noise and slight variations in 
the phase-setting translation velocity, it was more accu¬ 
rate in practice to determine the unknown phases 
by subtracting the discrete Fourier transforms of the out¬ 
put signals. 

The resulting 12 measured maps are nearly uni¬ 

tary within error: over all 12 circuits the maximum value 
of = — ljk\ was 0.050, with a mean of 

0.010, and on average Djk differed from zero by just 
1.6 standard deviations, as determined through Monte 
Carlo analysis using our uncertainties in r^®^® and 0^®^®. 
Those uncertainties were derived directly from the mea¬ 
sured variance in output power ratios and relative phase 
respectively taken over multiple trials. Comparison of 
the measured matrices U™®^® to U®^^® in Eq. (2) allows 
nearly direct determination of the net phase (j )n ; notably 
the values of 0^®^® are consistent with U®^^® to within er¬ 
ror. The splitting parameters On, n G {1,..., 4} can be de¬ 
termined from r™®^® using numerical optimisation. The 
results for and On for all 12 measured can¬ 

didate circuits are shown in Fig. 3. Notably the fidelity 
values shown were calculated directly from the measured 
matrices U™®^®, but agree to within error with the values 
calculated from the measured phase and reflectivity de¬ 
viations in the manner depicted in Fig. 2. The relatively 
higher variance in Os is due largely to the increased sen¬ 
sitivity of longer couplers to slight index mismatches, as 
explained in Sec. HB; BS3 has a coupling region almost 
six times the length of BSl and BS2 in order to achieve 
the required internal phase shift. 

For the best measured device, the mode and pro¬ 
cess fidelities determined were F^(U™®^®, = 

0.931T0.001 and Fp=0.680±0.005 respectively; the 
net internal phase found was 0w= — 0.346±0.013, 
and the splitting parameter deviations determined 
were A(9i=0.087±0.002, A(92=0.083±0.003, AO 3 = 
—0.065±0.003, and A6>4=0.337±0.002. As an illustration 
these parameter deviations are also depicted in Fig. 2 
along with their individual effects on the fidelities. For 
all 12 devices BS4 was erroneously fabricated with a re¬ 
flectivity near 60% rather than the ideal value of 90.8% 
. However we note that if A ^4 had been approximately 
the mean of the other splitting deviations achieved, with 
a value of 0.08, the process fidelity calculated according 
to those parameter errors would have been F^=0.882 and 
the mode fidelity F^ =0.984. 


B. Verification using two-photon interference. 

Full operation of these circuits as gates using sponta¬ 
neous parametric down-conversion (SPDC)—the current 
state-of-the-art in generating multiple single photons— 
requires a six-photon output state where two serve as 
triggers. This is to avoid heralding false positives due 
to the probabilistic nature of SPDC. Given the 2/27 gate 
success probability and our loss of at least 1.8 dB per 
waveguide, this would result in a success probability of 
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FIG. 3. Results from coherent circuit characterisation. Error bars are too small to see and are thus absent. Top: Mode and 
process fidelity of the measured circuit mappings with the ideal hcz circuit unitary, optimised over local external phases. The 
mean uncertainty in Fm and Fp are 0.002 and 0.006 respectively; these uncertainties were determined via Monte Carlo methods 
using the measured uncertainties in phases and moduli. Centre: Net undesired internal phase The displayed value is the 
mean of the four values determined from the four occurrences of in comparing for each circuit to Eq. (2). Eor all 

12 candidate circuits these differ by a maximum of 0.07. The variation between these four values dominated that between our 
many phase measurement trials for each circuit, and their standard deviation is thus taken to be our uncertainty; the mean 
resulting value over all 12 circuits is 0.0015. Bottom: Deviations AO from the ideal reflectivity parameters for the four BSs. Eor 
all values of AO the mean uncertainty—determined from repeated measurement trials—is 0.0026. In all 3 panels, the measured 
circuits are ordered by decreasing process fidelity. 


less than 0.015 per six-photon input. With current six- 
photon SPDC generation capabilities [39] we would thus 
expect a success rate of less than 3 mHz, necessitating 
prohibitively long integration times and low signal-to- 
noise for e.g. quantum process tomography, thereby lim¬ 
iting conclusions regarding actual gate fidelity. 

However, despite requiring four input photons for full 
operation, the hcz circuit relies only on fourth-order in¬ 
terference effects (in field), i.e. two-photon quantum in¬ 
terference; any higher-order interference effects between 
the four input photons can only result in error terms 
where the control and target qubits along with the two 
ancillary modes do not output exactly one photon each. 
In order both to confirm the ability of the circuits to sup¬ 
port high-visibility quantum interference and to verify 
the results of their coherent characterisation, we there¬ 
fore measured the visibility of two-photon quantum in¬ 
terference in the best-performing circuit for all possible 


input-output mode combinations. These measured visi¬ 
bilities are compared against both predictions from our 
classical characterisation and the ideal hcz circuit visi¬ 
bilities in Fig. 4. 

The apparatus for measuring the quantum interfer¬ 
ence effects is depicted schematically in Fig. 4(a). We 
measured two-photon quantum interference visibility for 
all ( 2 ) =36 combinations of two input and two out¬ 
put ports. The interference visibility V is calculated as 
V= {Cma^—Cmin) /C'max, whcrc C is the rate of coinci¬ 
dent photon detection events as a function of the tem¬ 
poral delay between the input photons, and Cmax and 
Cmin are calculated from a fit to the data as shown in 
Fig. 4(b). 

The measured visibilities are shown in Fig. 4(c), along 
with those predicted from as determined via co¬ 

herent characterisation, and the visibilities for an ideal 
hcz circuit. The mean of the residuals after subtract- 
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FIG. 4. (a) Setup for measuring two-photon interference. Degenerate photon pairs at 820 nm are created via spontaneous 
parametric down-conversion (SPDC) in a nonlinear /3-barium-borate (BBO) crystal pumped by a 410 nm frequency-doubled 
Ti:sapphire laser. Manual polarisation controllers (POL) enable alignment of SPDC polarisation with the axes of polarisation- 
maintaining fibers coupled to the test device. We detect photons in coincidence using avalanche photo diodes (APD). (b) 
Representative observed non-classical interference patterns, showing anti-coalescent and coalescent interference for two different 
output mode combinations as well as fits to the data with Gaussian and sine components. The error bars shown stem from 
Poissonian counting statistics, (c) Two photon interference results for our best hcz circuit. We compare predictions (left) 
from the coherently-characterised circuit against measured two-photon interference visibilities (centre). The right panel shows 
visibilities for the ideal circuit most of the difference between this panel and the other two is due to the deviation in 

our best circuit from the ideal reflectivity for BS4 and unwanted net phase as shown in Fig. 3. For the residuals after 
subtracting the measured visibilities from the predictions the mean and standard deviation are -0.002 and 0.061, while after 
subtracting the measured values from the ideal values they are -0.022 and 0.282 respectively. 


ing the measured visibilities from the predictions is only 
-0.002 with a standard deviation of 0.061. Perhaps a bet¬ 
ter comparison is achieved by numerically calculating the 
unitary which would yield the minimum root-mean- 
square difference from the measured visibilities; this uni¬ 
tary has mode fidelities of = 0.983 and 

= 0.931 with the measured circuit and 
the ideal hcz respectively. The small differences between 
predicted and measured visibilities can be attributed 
largely to three factors: polarisation non-degeneracy be¬ 
tween the interfering photons in the FLDW circuit; the 
slightly differing spectra of the SPDC photons and the 
laser diode used for the coherent characterisation; and 
the effects of higher-order SPDC terms. 


IV. CONCLUSIONS 

Along with further improvements in photon sources 
and detection, heralding will be required to concate¬ 
nate multiple entangling LOQC gates and thus enable 
more complex quantum computations and simulations. 
We have demonstrated that integrated waveguide arrays, 
particularly using femtosecond laser-writing, are capable 
of generating the required multimode interference circuits 
with both high fidelity and excellent quantum interfer¬ 
ence, and allow simple implementation of mode crossover 
elements and tt phase shifts. However further careful en¬ 
gineering will be required to precisely achieve the desired 
beamsplitter ratios and to avoid undesired phase accu¬ 
mulation. This study has outlined the challenges both 
experimentally and theoretically in achieving waveguide 
circuits with high operational fidelities. 

The quantum process fidelity of candidate circuits can 
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be calculated from known fabrication tolerances or clas¬ 
sical characterisation results using the Jamiolkowski iso¬ 
morphism, and this metric has proven to be more sen¬ 
sitive and useful than mode fidelity for assessing such 
circuits. However any circuit imbalance will lead to error 
terms outside the computational subspace wherein two 
photons exit in either the control or target mode. The 
precise effects of such coherent error terms when multiple 
gates are concatenated, as well as possibilities for their 
correction or mitigation, could be a fruitful avenue for 
future investigation. 
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